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Abstract. - We consider a simple model for examining the effects of quenched disorder on drag 
consisting of particles interacting via a Yukawa potential that are placed in two coupled one- 
dimensional channels. The particles in one channel are driven and experience a drag from the 
undriven particles in the second channel. In the absence of pinning, for a finite driving force 
there is no pinned phase; instead, there are two dynamical regimes of completely coupled or 
locked flow and partially coupled flow. When pinning is added to one or both channels, we find 
that a remarkably rich variety of dynamical phases and drag effects arise that can be clearly 
identified by features in the velocity force curves. The presence of quenched disorder in only the 
undriven channel can induce a pinned phase in both channels. Above the depinning transition, 
the drag on the driven particles decreases with increasing pinning strength, and for high enough 
pinning strength, the particles in the undriven channel reach a reentrant pinned phase which 
produces a complete decoupling of the channels. We map out the dynamic phase diagrams as 
a function of pinning strength and the density of pinning in each channel. Our results may be 
relevant for understanding drag coupling in ID Wigner crystal phases, and the effects we observe 
could also be explored using colloids in coupled channels produced with optical arrays, vortices in 
nanostructured superconductors, or other layered systems where drag effects arise. 
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There are many examples of one- and two-dimensional 
('ID and 2D) coupled bilayers or coupled channels of in- 
teracting particles, including vortices in superconducting 
bilayers colloidal systems dusty plasmas (T], 

and Wigner crystals [8HTU] . In many of these systems it is 
possible to apply an external drive to one of the layers and 
ipeasure the resulting response of the other layer as well as 
the drag effect produced by the particles in the undriven 
l^,yer on those in the driven layer. For instance, this type of 
measurement has been performed for the transformer ge- 
ometry in two-layer superconducting vortex systems [Tl[2] . 
If the vortices in each layer are completely coupled, the 
measured response is the same in both layers. If instead 
the vortices are only partially coupled between layers, the 
response is reduced in the undriven layer compared to 
the driven layer. A similar coupling-decoupling transition 
is also predicted for coupled wires containing ID Wigner 
crystals [5], and certain drag effects in ID wires have been 
interpreted as resulting from the formation of ID Wigner 
crystal states [5] . 



Many of these systems can contain some form of 
quenched disorder which could produce pinning effects 
however, very little is known about how pinning al- 
ters the drag or transport properties in layered geometries. 
The quenched disorder could be strong in one channel and 
weak in another, or it could be of equal strength in both 
channels, resulting in different types of dynamic phases. 
Particles driven over random quenched disorder in single 
layer systems are known to exhibit a variety of dynamical 
phases associated with distinct transport signatures, as 
shown in simulations jllH14j and experiments [15| . These 
studies suggest that distinct dynamics could also arise in 
a two- layer drag system when pinning is present. We note 
that there have been studies of two-layer driven colloidal 
systems; however, the effect of quenched disorder was not 
considered [4]. 

In this work we propose a simple model of two coupled 
ID channels of particles with repulsive Yukawa interac- 
tions where there a drive channel and a drag channel, and 
where pinning is added to one or both channels. Despite 
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Fig. 1: A schematic of the system of two coupled ID channels 
separated by a distance d. Each channel contains A*' particles 
(black dots) which interact via a Yukawa potential with parti- 
cles in the same and the adjacent channel. An external drive 
Fr> is applied to the particles in the drive channel, resulting in 
an average particle velocity of Vdrive in the drive channel. The 
velocity response in the drag channel is Vdrag- Quenched dis- 
order is introduced in the form of localized pinning sites (open 
circles) of maximum force Fp which are spaced randomly in the 
channel without overlapping. The pinning can be placed only 
in the drag channel (as shown), only in the drive channel, or 
in both channels. 



the apparent simplicity of this system, we find that a re- 
markably rich variety of distinct dynamical phases are pos- 
sible which produce pronounced changes in the transport 
and drag effects. Although our model treats ID channels, 
we expect that many of the same effects should be generic 
to coupled layer systems with quenched disorder. The spe- 
cific system we consider could be realized experimentally 
using colloidal particles in channel geometries |16j where 
a driving force is applied to one channel via optical means 
or an electric field. The effects we observe could also be 
studied with coupled ID Wigner crystals [9ll 11 II 171118) or in 
superconductors with nanofabricated channel geometries 
[inj when the current is applied to only one channel and 
the response is measured in both channels. In our system, 
we find that without pinning there is a finite drive transi- 
tion from a locked regime to a partially decoupled regime 
in which the response of the drag channel decreases with 
increasing drive. Addition of pinning to the drag chan- 
nel can induce a pinned phase for both channels. As the 
drive increases, both channels depin into a locked or par- 
tially locked phase, followed at high drives by a reentrant 
pinning of the drag channel when the dynamic coupling of 
the two channels becomes weak enough. When the driving 
force is fixed at a low value, increasing the strength of the 
pinning in the drag channel increases the effective drag on 
the particles in the driven channel; for higher fixed driving 
force, however, increasing the pinning strength reduces the 
drag on the driven channel and causes the channels to de- 
couple when the drag channel becomes reentrantly pinned. 
We also observe strong fluctuations in both channels at 
the locked to partially locked transition for intermediate 
strengths of quenched disorder. We map out the evolution 
of the dynamical phases for a wide variety of parameters 
including pinning in both channels, pinning in only one of 
the channels, and for varied pinning densities. 

Simulation- In fig. 1 we show a schematic of our system 
which consists of two ID channels separated by a distance 
d with periodic boundary conditions in the a;-direction. 
Each channel contains N particles which interact via a re- 



pulsive Yukawa potential with other particles in the same 
channel and with particles in adjacent channels. Particles 
in the upper or drive channel are subjected to an applied 
external force Fd. A single particle i located at position 
R; undergoes motion obtained by integrating the over- 
damped equation of motion: 



V- 



'dt 



(1) 



Here ry 
to f] 
tion 



is the damping constant which is set 
= 1.0. The particle-particle interac- 



forcc is FPP = Ej^I, ) where 
V{Rij) = {Eq / Rij) exp{— KRij)'R.ij , Rij = |Ri — Rj|, 
Rij — (Ri — Rj)/i?ij, and Eq — Z*'^ /Aireeoao. Here e is 
the solvent dielectric constant, Z* is the effective charge, 
and 1/k is the screening length. In each channel, the 
average particle spacing is a and we take d/a = 2/3. 
We also take 1/k = 2d to ensure that coupling between 
the two channels is possible. We note that we have 
also considered other values of d, a, and k and find the 
same qualitative features, indicating that our results 
should be generic for this class of system. The pinning 
force Ff arises from Np attractive parabolic potentials 
with a maximum force of Fp and a radius of Rp ~ 0.5, 

= EkZi EoFpiRP,/Rp)QiRp - R.,k)-RP,. Here 9 is 
the Heavisidc step function, Rf, is the location of pinning 
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this work we consider the limit in which the number of 
pinning sites is smaller than or similar to the number 
of particles in the system; however, we find the same 
qualitative features when we change the pinning density 
or strength. The external force Fd is applied only to the 
particles in the driven channel. We increase Fo from zero 
in small increments of SFd = 0.001, spending a waiting 
time of 10** simulation time steps at each increment. We 
measure the average velocity of the particles in the drive 
and drag channels, Vdrag = ^d'rig SilT"" '^^i/^* ^^'^ 
Vdrive = N'^l^^Y^fjl"^ dxi/dt, where Ndrag = N is the 
number of particles in the drag channel and Ndrive = N 
is the number of particles in the drive channel. This 
simulation method, employing overdamped dynamics 
with random pinning, has been used previously to explore 
the behavior of colloids driven over random and periodic 
substrates as well as two layer systems. A similar method 
has been used to study sliding Wigner crystals and 
vortices in type-II superconductors. 

In fig. 2(a) we show the velocity force curves Vdrive and 
Vdrag vcrsus Fu foT the case where there is no pinning 
in either channel. At low _Fd, Vdrive — Vdrag, indicat- 
ing that the particles are moving together in the coupled 
region (C). At Fu = 2.0, the average particle velocity 
in each channel is just under 1.0. In the absence of the 
drag channel, we would have Vdrive = 2.0 at Fd = 2.0 
since V = Fu/rj under our overdamped dynamics. When 
the drive channel is locked with the drag channel, the 
same driving force must be used to transport twice as 
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Fig. 2: The average particle velocity in the drive channel Vdrive 
and the drag channel Vdrag vs applied force Fd for a sample 
with pinning in the drag channel at Np/N = 0.083. (a) The 
pin-free case Fp — 0.0 has two regions: an initial coupled re- 
gion (C) where motion in both channels is completely locked, 
and a partially coupled (PC) region where Vdrag monotoni- 
cally decreases while Vdrive increases with increasing Fd- (b) 
For Fp = 1.0 there is an initial pinned region (P) in both chan- 
nels in addition to the C and PC regions, (c) At Fp = 2.5, 
an intermittent region (/C) containing a series of asymmetric 
jumps in Vdrag appears between the C and the PC regions. 
At Fd = 5.66 there is a reentrant pinning of the drag chan- 
nel above which Vdrag ~ 0.0, producing an uncoupled region 
(UC). (d) At Fp = 4.0, only the P and UC regions appear and 
Vdrag = for aU Fd ■ 



many particles, reducing the particle velocity by half. At 
Fd = 2.124 there is a transition to a partially coupled 
(PC) state where the velocity curves branch apart. Vdrive 
continues to increase with increasing Fjy with a slope 
higher than the slope of Vdrive below Fd = 2.124, while 
Vdrag dccreases with increasing Fjj. Since there is no pin- 
ning, in the PC region Vdrag gets smaller and smaller with 
increasing F^ but never reaches zero. The general shape 
of the velocity force curves is very similar to the response 
observed for the superconducting transformer geometry 
where vortices in a pair of adjacent superconducting lay- 
ers can couple and decouple depending on the applied cur- 
rent [HH]. In the superconducting system, the voltage 
response is proportional to the particle velocity and the 
applied current is proportional to Fr). Additionally, the 
velocity-force curves in fig. 2(a) are in agreement with the 
predicted coupling-decoupling transition for disorder-free 
coupled ID Wigner crystal systems [5]. 

In fig. 2(b) we plot Vdrive and Vdrag versus Fjj for 
a system containing pinning in the drag channel with 
Np/N ~ 0.083 and Fp = 1.0. The curves are similar to 
the pin free case, but a pinned region (P) now appears for 
Fjj < 0.1 when the quenched disorder in the drag channel 
causes the particles in both channels to stop moving. The 
particles in the drag channel are directly pinned by the 



pinning sites, but the particles in the driven channel are 
only indirectly pinned by their interactions with the par- 
ticles in the drag channel, which produce an effective peri- 
odic pinning potential for the driven particles. The driven 
particles can either simply move over this periodic poten- 
tial, in which case we find a decoupled region {UC) where 
the driven particles are moving and the drag particles are 
pinned, or the driven particles can drag the periodic poten- 
tial along with them, depinning the drag particles and pro- 
ducing a moving coupled region. As Fjj increases, there 
is a transition from the coupled to partially coupled flow 
at the same Fo = 2.124 as in the pin- free system. Fig- 
ure 2(c) shows the velocity- force curves for a sample with 
Fp = 2.5, where the transition out of region C falls at a 
lower Pd = 1.26 and is followed by a region containing a 
series of jumps in both Vdrive and Vdrag- In this intermit- 
tently coupled (/C) region, each downward jump in Vdrag 
is preceded by a range of F^ over which Vdrag increases 
linearly with increasing Fjj. The jumps in Vdrag arc ac- 
companied by similar sharp upward jumps in the value of 
Vdrive- As Fd increases, the jumps decrease in size un- 
til the system crosses over to the partially coupled (PC) 
region where Vdrag decreases linearly and smoothly with 
increasing Fd- At high Fd, the uncoupled channel {UC) 
region appears when the particles in the drag channel be- 
come pinned again for Fd > 5.66 but the particles in the 
driven channel continue to flow. The reentrant pinning 
of the drag channels results when the effective dynami- 
cal coupling between the two channels weakens for higher 
Fd, as indicated by the decrease in Vdrag with increas- 
ing Fd in the PC region of the pin-free sample shown 
in fig. 2(a). When the dynamical coupling between the 
channels drops below a critical value, the drag force ex- 
perienced by the particles in the drag channel falls below 
the depinning threshold for the pinning sites in the drag 
channel, resulting in the repinning transition. We expect 
this reentrant pinning effect to be a generic feature in drag 
systems containing quenched disorder, and its signature is 
a sudden decrease of the drag channel response to zero. 
In fig. 2(d) for Fp = 4.0, the drag channel never depins 
over the range of Fd shown, and there is a single transi- 
tion from region P where both channels are pinned to the 
decoupled region UC- 

In fig. 3(a) we plot the dynamic phase diagram Fd ver- 
sus Fp for the system in fig. 2 containing pinning only 
in the drag channel. The transition between the PC and 
U C regions drops to lower values of Fd with increasing Fp 
since a lower drive is required to reduce the effective drag 
force below the depinning threshold of the drag channel 
at higher Pp. For Fp > 3.88, only the P and UC regions 
appear. The extent of region C decreases with increasing 
Fp until region C disappears for Fp > 3.45. The IC region 
vanishes at higher values of Fd and reaches its maximum 
extent near Fp = 2.95. In fig. 3(b) we plot Vdrive versus 
Fp for the same system at fixed Fd = 0.2. Here Vdrive 
decreases with increasing Fp until the system enters the 
pinned region at Fp = 2.0 and Vdrive drops to zero. The 



p-3 



C. Bairnsfather, C.J.O. Reichhardt and C. Reichhardt 




Fig. 3: (a) The dynamic phase diagram of Fd vs Fp for the sys- 
tem in fig. 2 for pinning in only the drag channel. Each region 
is marked: pinned (P), coupled (C), intermittently coupled 
(IC), partially coupled (PC), and uncoupled (UC). (b) Vdrive 
versus Fp from the system in (a) at Fd = 0.2. The increas- 
ing effective damping 77* = Fd /Vdrive produces a monotonic 
decrease in Vdrive with increasing Fp until the particles in the 
driven channel repin and the velocity drops to zero in region 
P- (c) Vdrive VS Fd for the system in (a) at Fd ~ 3.0. Here 
Vdrive increases with increasing Fp until the onset of the UC 
phase, above which Vdrive saturates to a constant value. 



effective damping of the driven particles, rj* ~ Fd /Vdrive, 
increases with increasing Fp. In fig. 3(c) we show Vdrive 
versus Fp at fixed Fjj = 3.0, where Vdrive monotonically 
increases until the system enters the UC phase, and Vdrive 
saturates to a constant value. These results show that pin- 
ning in the drag channel can either increase or decrease 
the effective drag on the driven channel depending on the 
strength of the pinning and the amplitude of the dc drive. 

We next consider a system containing pinning in both 
channels. Figure 4(a) shows the dynamic phase diagram 
of F]j versus Fp for a sample with a pinning density of 
Np/N — 0.083 in each channel. The overall shape of the 
phase diagram resembles the phase diagram in fig. 3(a) 
obtained with pinning in only the drag channel. The 
pinned phase P grows in size with increasing Fp much 
more rapidly when there is pinning in both channels than 
when there is pinning in only the drag channel. In fig. 4(c) 
we plot Vdrive and Vdrag versus Fjj at Fp = 4.0 for the sam- 
ple from fig. 4(a) with pinning in both channels. When 
Vdrag drops to zcro at Fjj = 2.22 at the transition to re- 
gion UC, there is a jump up in the value of Vdrive- In 
fig. 4(b) we show the phase diagram of Fjj versus Fp for a 
sample containing pinning in only the drive channel with 
Np/N = 1.33. The UC phase no longer appears since the 
particles in the drag channel can only reentrantly repin 
when there is pinning in the drag channel. For Fp > 5.5, 
the system crosses directly from region P to region PC 
without passing through the C region. Figure 4(d) il- 
lustrates Vdrive and Vdrag vcrsus F]j for this system at 
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Fig. 4: (a) The dynamic phase diagram of Fd vs Fp for 
a system with pinning in both channels. Each channel has 
Np/N — 0.083. The region labels are the same as in fig. 3(a). 
(b) The dynamic phase diagram of Fd vs Fp for a system with 
pinning in only the drive channel, with Np/N = 1.33. Here 
the uncoupled channel region is absent and there can be a di- 
rect transition from region P to region PC. (c) Vdrive (upper 
curve) and Vdrag (lower curve) vs Fd for the system in (a) at 
Fp = 4.0. Inset: Vdrive vs Fp for the system in (a) at Fd = 3.0. 
(d) Vdrive (upper curve) and Vdrag (lower curve) vs Fd for the 
system in (b) at Fp = 7.5. Inset: Vdrive vs Fp for the system 
in (b) at Fd = 3.0. 



Fp ~ 7.5. There is a transition directly from region P 
to region PC at Fd = 4.8, and above this drive Vdrag 
gradually decreases with increasing Fd ■ 

In the inset of fig. 4(c) we plot Vdrive versus Fp at Fd = 
3.0 for the system in fig. 4(a) where there is pinning in both 
channels. For Fp < 3.5, Vdrag increases with increasing 
Fp, just as shown in fig. 3(c) for the sample with pinning 
in only the drag channel. For Fp > 3.5, when there is 
pinning in both channels Vdrag decreases with increasing 
Fp for Fp > 3.5 rather than saturating. The increase in 
Vdrag with increasing Fp for Fp < 3.5 results when the 
pinning in the drag channel reduces the coupling between 
the driven and drag particles. For Fp > 3.5 at Fd = 3.0, 
the particles in the drag channel are always pinned. If 
there is pinning in only the drag channel, Vdrive remains 
constant for Fp > 3.5 as shown in fig. 3(c); however, if 
pinning is also present in the drive channel, it produces 
an increasing drag on the driven particles, causing Vdrive 
to decrease with increasing J^p at Fp > 3.5 for the sample 
with pinning in both channels. In the inset of fig. 4(d) we 
illustrate this effect more clearly by plotting Vdrive vs Fp 
for a system with pinning in only the drive channel. Here 
Vdrive monotonically decreases with increasing Fp until it 
reaches zero once all the drive particles become pinned. 

It is also possible to vary d/a, where d is the distance 
between the channels and a is the average particle spac- 
ing within the channels. We consider the effect of varying 
particle density in samples with fixed d. Increasing the 
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Fig. 5: (a) The dynamic phase diagram of Fd vs d/a for a 
sample with pinning in both channels at Fp = 2.0 and with 
Np/N = 0.25 at d/a — 0.67. (b) The dynamic phase diagram 
of Fd vs Np/N for a sample with pinning in both channels at 
d/a = 0.67 and Fp = 1.0. 



particle density reduces the effective coupling between the 
channels since the periodic potential experienced by the 
particles in one channel produced by the particles in the 
other channel becomes smoother as the particle density 
increases. In fig. 5(a) we plot the dynamic phase dia- 
gram for Fd versus d/a in a sample with pinning in both 
channels for fixed Fp = 2.0 and with Np/N — 0.25 at 
d/a = 0.67. At low d/a, there is a large pinned region 
since the particles within each channel are far from each 
other, strongly enhancing the effectiveness of an individual 
pinning site. For d/a < 0.22, the system passes directly 
from region P to region UC. For large d/a, each particle 
interacts strongly with other particles in the same chan- 
nel but the coupling to particles in the other channel is 
greatly reduced, so the system again crosses directly from 
the pinned to the UC region. For intermediate values of 
d/a, the three other regions C, IC, and PC appear in win- 
dows between the P and UC regions. The total width of 
these intermediate windows reaches a peak at d/a — 0.66. 
In fig. 5(b) we show the dynamic phase diagram of F]j ver- 
sus the pinning density Np/N for a system with pinning 
in only the drag channel at Fp = 1.0 and d/a = 0.67. The 
transition from region PC to region UC drops to lower Fu 
with increasing Np/N since higher pinning density causes 
the repinning threshold in the drag channel to decrease. 



The transition between region C and region PC also shifts 
to lower F]j with increasing Np/N, and there is a narrow 
window (not labeled on the figure) where region IC ap- 
pears between regions C and PC. There is also a small 
plateau in all of the transitions near Np/N = 1.0, which is 
a weak comnicnsuration effect. For periodically arranged 
pinning sites, we expect much stronger commensuration 
effects to occur; these effects will be considered elsewhere. 

In summary, wc introduce a simple model to study the 
effects of pinning on drag in coupled one-dimensional chan- 
nels of particles interacting via a Yukawa potential where 
a drive is applied to only one channel. This system ex- 
hibits a rich variety of dynamical phases which are asso- 
ciated with distinct features in the velocity-force curves. 
In the absence of pinning, at low drives both channels 
are coupled and exhibit identical particle velocities. As 
the drive is increased, there is a transition to a partially 
coupled phase where the velocity in the drag channel de- 
ceases with increasing drive while the velocity in the driven 
channel increases with increasing drive. Placing pinning 
in only the drag channel produces a region in which both 
channels are pinned as well as a reentrant pinning of the 
particles in the drag channel at higher drives when the 
coupling between the channels is dynamically weakened. 
When strong pinning is placed either in the drag channel 
or in both channels, we observe only a pinned region and 
a decoupled region. For intermediate pinning strength, we 
find a strongly fluctuating region where a series of asym- 
metric jumps appear in the velocity-force curves. Increas- 
ing the strength of the pinning in the drag channel gen- 
erally reduces the effective drag on the particles in the 
driven channel at higher drives since the pinned particles 
in the drag channel are unable to absorb momentum from 
the particles in the drive channel. We map the dynamic 
phase diagrams for samples with pinning in both channels, 
pinning in only the drive channel, and pinning in only the 
drag channel. The smallest number of phases appears for 
pinning in only the drive channel, while samples with pin- 
ning in both channels exhibit multiple dynamical phase 
transitions even when the particle density or pinning den- 
sity is altered, showing that the effects we observe are ro- 
bust for a wide range of parameters. Our results could be 
important for understanding drag effects in systems where 
ID Wigner crystal states can occur and for understanding 
the general effects of quenched disorder in other drag and 
transformer geometries. Additionally, our system could 
be directly realized experimentally for colloidal or dusty 
plasma systems in coupled ID channels. Extensions of this 
model that could be explored include the effects of a peri- 
odic pinning array at and away from commensurability, as 
well as the effects of disorder on the coupling of particles 
confined in adjacent two-dimensional planes. 

^ ^ ^ 
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